arXiv:1505.07286vl [astro-ph.GA] 27 May 2015 


A Herschel/YUFl Legacy Survey of HF and H2O in the Galaxy: 
Probing Diffuse Molecular Cloud Chemistry^ 

P. Sonnentrucker^ 

Space Telescope Science Institute, Baltimore, MD 21218/European Space Agency 

sonnentrQstsci.edu 

M. Wolfire^ 

University of Maryland, College Park, MD 20742 
D. A. Neufeld^^ 

Johns Hopkins University, Baltimore, MD 21218 

N. Flagey"^ 

Institute for Astronomy, Hilo, HI 96720-2700 
M. Gerin®’® 

Sorbonne Universites, Universite Pierre et Marie Curie, Paris 6, CNRS, Observatoire de 

Paris, UMR 8112, LERMA, Paris, France 

LERMA, Observatoire de Paris, PSL Research University, CNRS, UMR 8112, F-75014, 

Paris, France 

P. Goldsmith^ 

Jet Propulsion Laboratory, 4800 Oak Grove Drive, Pasadena CA 91109-8099 

D. Lis5’6,8 

Sorbonne Universites, Universite Pierre et Marie Curie, Paris 6, CNRS, Observatoire de 

Paris, UMR 8112, LERMA, Paris, France 

LERMA, Observatoire de Paris, PSL Research University, CNRS, UMR 8112, F-75014, 

Paris, France 


2 


and 

R. Monje® 

California Institute of Technology, Cahill Center for Astronomy and Astrophysics 301-17, 

CA 91125, Pasadena, USA 

Received_; accepted_ 



- 3 - 

ABSTRACT 

We combine Herschel observations of a total of 12 sources to construct the 
most uniform survey of HF and H 2 O in our Galactic disk. Both molecules are 
detected in absorption along all sight lines. The high spectral resolution of the 
Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the 
HF and H 2 O distributions in 47 diffuse cloud components sampling the disk. 
We hnd that the HF and H 2 O velocity distributions follow each other almost 
perfectly and establish that HF and H 2 O probe the same gas-phase volume. Our 
observations corroborate theoretical predictions that HF is a sensitive tracer of 
H 2 in diffuse clouds, down to molecular fractions of only a few percent. Using 
HF to trace H 2 in our sample, we hnd that the iV(H20)-to-iV(HF) ratio shows a 
narrow distribution with a median value of 1.51. Our results further suggest that 
H 2 O might be used as a tracer of H 2 -within a factor 2.5- in the diffuse interstellar 
medium. We show that the measured factor of ~2.5 variation around the median 
is driven by true local variations in the H 2 O abundance relative to H 2 throughout 
the disk. The latter variability allows us to test our theoretical understanding 
of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that 
both gas-phase and grain-surface chemistry are required to reproduce our H 2 O 
observations. This survey thus conhrms that grain surface reactions can play a 
signihcant role in the chemistry occurring in the diffuse interstellar medium (uh 
< 1000 cm“^). 

Subject headings: ISM: lines and bands — ISM: molecules — ISM: clouds — ISM: 
cosmic rays — ISM: abundances — Physical Data and Processes: astrochemistry 
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1. Introduction 

Diffuse molecular clouds are regions in which atomic hydrogen is progressively 
converted to molecular hydrogen (H 2 ), and neutral carbon and carbon monoxide become 
the dominant forms of carbon as the total visual extinction increases from 0.1 (the onset 
of H 2 formation) to about 2 magnitudes (see, Snow & McCall 2006, for a review). Diffuse 
molecular clouds are considered clouds in transition from diffuse mainly atomic gas to fully 
molecular gas, hence an important first step in star formation. As a result, they play a 
crucial role in the lifecycle of the interstellar medium making their study critical to advance 
our understanding of how molecular clouds form from the diffuse interstellar medium (ISM). 

Because of their relatively low densities (uh <1000 cm“^) and their low shielding 
from UV radiation compared to dense clouds, diffuse molecular clouds were expected 
to be mostly devoid of molecules. However, the last four decades of UV/optical/radio 
observations, from space and from the ground, have demonstrated that diffuse molecular 
clouds have a surprisingly rich and still largely unexplained chemistry (e.g.. Snow & McCall 
2006, Sonnentrucker et al. 2007; Sheffer et al. 2008; Liszt et al. 2007; Neufeld et ah 2012). 
Comparisons of UV/optical molecular absorption line studies with high-resolution studies 
using sub-mm data have provided complementary information on the physics and chemistry 
of the absorbing gas, over a large range of opacities and at great distances in the Galactic 
disk. As a result, these clouds constitute “in-situ” laboratories in which we can study a 
variety of physical and chemical processes of broad applicability in astrophysics. 

Small (diatomic and triatomic) hydrides are important tracers of the diffuse interstellar 
medium physics and chemistry. Possessing small moments of inertia, these hydrides 

^Herschelis an ESA space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from NASA. 
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have rotational transitions at THz frequencies that are difficult or impossible to observe 
using ground-based observatories. The Herschel/Y{\¥\ Key Program PRISMAS (Probing 
Interstellar Molecules with Absorption line Studies: PI, M. Gerin) was aimed at surveying 
key hydrides within the Galaxy. Up to 22 small molecular species were specifically targeted 
with PRISMAS in order to probe the chemistry of carbon (e.g., Mookerjea et al. 2010, 
Gerin et al. 2010, Godard et al. 2012; 2014), nitrogen (e.g., Persson et al. 2010; Persson 
et al. 2014), oxygen (e.g., Neufeld et al. 2010; Flagey et al. 2013; Indriolo et al. 2015), 
chlorine (e.g., Lis et al. 2010; De Luca et al. 2012; Neufeld et al. 2012, Monje et al. 2013) 
and fluorine (e.g., Neufeld et al. 2010, Sonnentrucker et al. 2010; this paper) in diffuse 
molecular clouds and to constrain the physical processes at play in the Galactic diffuse ISM. 

In this work, we focus on two particular hydrides, hydrogen fluoride (HP) and water 
(H 2 O) with the aim of further probing the chemistries of fluorine and oxygen-bearing 
molecules and of determining to what extent both species can be used as diagnostics of 
the physical processes at play in the diffuse ISM. The combination of the full PRISMAS 
sample, with some WISH data (Water In Star-forming regions with Herschel; PL E. F. van 
Dishoeck) and Herschel Gycle 1 data allows us to report on the largest set of HF and H 2 O 
column density measurements in the Galactic disk to-date. Sections 2 & 3 describe our 
set of observations, as well as the reduction and analysis techniques we employed. Section 
4 summarizes our results. In Sections 5 & 6 we compare the HF and H 2 O distributions 
measured with our survey to the most recent 2-sided PDR model predictions (Neufeld & 
Wolhre 2009, Hollenbach et al. 2012) for the range of physical conditions most relevant to 
diffuse molecular clouds. In particular, we validate the role of HF as a tracer of H 2 and 
add weight to the suggestion that H 2 O can be used as a tracer of H 2 in diffuse clouds as 
well (e.g., Flagey et al. 2013). We conhrm the co-spatial distribution of HF and H 2 O in 
the Galactic diffuse ISM, which in turn allows us to derive the H 2 O abundance relative 
to HF throughout the Galactic disk. In Section 7 we discuss how the measured variations 



in the H 2 O abundance relative to HF can be used to test our understanding of the H 2 O 
chemistry. We assert the role that grain surface chemistry plays in the production of H 2 O 
in the diffuse interstellar medium. Section 8 summarizes our general conclusions. 


2. Observations 

The names and coordinates of the 12 targets comprising our sample are summarized 
in Table [H All targets are well-studied star-forming regions, known to produce strong 
far-infrared background continuum emission, a prerequisite for detecting absorption by 
foreground interstellar clouds. The sight lines to each of these sources are also known 
to intersect multiple spiral arms, hence, allowing us to probe a variety of local physical 
conditions throughout the Galactic disk (e.g., Godard et ah 2010; Lang et ah 2010). 

The sources G—0.02—0.07 (-1-50 km s~^ cloud around SgrA), G—0.13—0.08 (-1-20 km 
s-i cloud around SgrA), W28A, W31G, W33A, G34.3+0.1, W49N, W51, W3 IRS5 and 
DR21(OH) were observed through Guaranteed Time Key Program PRISMAS. W3 IRS5 
was also observed in Guaranteed time Key Program WISH. W3(OH) and G29.96—0.02 were 
observed as part of the Herschel Open Time Gycle 1 campaign (OTl, PI: D.A. Neufeld) 
along with W51 that was reobserved to signihcantly increase our signal-to-noise (S/N) ratio 
in the continuum and push our detection limits to lower column density clouds compared 
to those measured from data previously obtained with the PRISMAS program. 

We observed the ground-state rotational line of HF (rarest = 1232.476 GHz; Nolt et 
al. 1987) in the upper sideband of HIFI band 5a receiver and the ground-state line of 
para-H20 (P-H 2 O; v^est = 1113.343 GHz) in the lower sideband of HIFI band 5a receiver. 
We used multiple Local Oscillator (LO) settings in order to securely identify the HF and 
P-H 2 O absorption lines along all 12 sight lines. The proximity of the HF and P-H 2 O lines in 
frequency and their detection within the same receiver band ensure that the absorption lines 
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are observed with a very similar telescope performance and that the data are calibrated in 
a similar fashion. 

We nsed the Dual Beam Switch (DBS) mode which, combined with the Wide Band 
Spectrometer (WBS), allows for a spectral resolution of about 1.1 MHz (0.3 km s“^ at the 
HF frequency). The DBS mode uses two reference OFF-beam positions located 3' on either 
side of the source position, along an East-West axis. Because the Galactic Center is a very 
complex and crowded region, we checked for emission or absorption (contamination) in the 
OFF-beam position using Herschel observations performed by the HEXGAL project (R. 
Glisten and M. Requena-Torres, private communication). There is no contamination at 
the OFF-beam position for HF; only minor contamination occurs for P-H 2 O in the velocity 
range [— 40, -|-10] km s“h As a result, the absorption components detected in the Galactic 
disk in this particular velocity range are not included in the analysis described in Section 3 
for the two Galactic Center sources G—0.02—0.07 and G—0.13—0.08. Contamination in the 
OFF-beam position is not a concern for the remaining more compact sources in our sample. 


3. Data Reduction and Analysis 

The data were processed to Level 2 with the standard HIFI pipeline in the Herschel 
Interactive Processing Environment (HIPE) version 9.1 (Ott 2010), thereby producing fully 
calibrated spectra for both polarization modes at each LO setting. Further inspection of 
the Level 2 data showed that the signals obtained in each of the two polarization modes for 
a given LO setting were in excellent agreement, as were the HF and P-H 2 O spectra obtained 
at each LO setting. Occasionally, emission from molecules other than HF or P-H 2 O appear 
in the sideband containing the HF or P-H 2 O absorption (see Fig[T]). For 6 of the 12 sources 
presented here, Flagey et al. (2013) performed a detailed study of the distribution of water 
in the ground and excited states accessible to Herschel/HIFI. While contamination needed 


to be taken into account for some of the excited states of water, their study showed that 
the ground-state of P-H 2 O, of interest here, was mostly free of this effect. For the remaining 
6 sources in our sample we compared the absorption profiles of the ground states of HF 
and P-H 2 O across LO settings and for both polarizations. Our comparisons revealed no 
signihcant contamination with interloping features for both molecules. As a result, for each 
target, we generated average spectra for HF and P-H 2 O that consist of the weighted sum 
of up to 6 spectral observations (up to three LO settings with two polarizations each), for 
which each observation is weighted in inverse proportion to the square of its root mean 
square (rms) noise. The double sideband continuum antenna temperatures T^(cont) 
and respective rms noise for HF and P-H 2 O derived from these weighted average spectra 
are reported in Table [21 which also lists the total on-source exposure times (texp), the 
observation dates, and the observation IDs (Obs IDs) of the LO settings obtained for HF 
and P-H 2 O. 


3.1. Background Continuum Emission Treatment 

HIFI employs double sideband receivers and for a sideband gain ratio equal to unity, 
the saturated absorption of radiation at a given frequency for a transition with excitation 
temperature much less than TA(cont) will reduce the measured antenna temperature (Ta) 
to one-half the apparent continuum antenna temperature TA(cont). 

The top two panels of Figures 1 through 6 display the double sideband weighted average 
spectra for HF and P-H 2 O, respectively, versus Doppler velocity in the Local Standard of 
Rest frame (Vlsr) for each of the 12 sight lines probed in our survey. For most sight lines, 
HF is detected only in absorption, and the background continuum emission from the source 
itself is well modeled with a constant (red line in each panel). The flux normalized with 
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respect to the continuum flux in a single sideband can then be expressed as, 

Tyi/TA(cont) = (1/2) [1 + exp(-r)] 

assuming that the sideband gain ratio is equal to unity. The latter assumption was 
investigated by Flagey et ah (2013) who found that the HIFI sideband gain ratio in band 
5a, the band of interest here, was very stable and indeed consistent with unity. 

Toward W31C, W3 IRS5, W28A and W49N, one can see that HF emission arising 
from gas local to the star-forming region blended itself with HF absorption arising from 
foreground clouds with projected velocities similar to those of the emitting gas. Such 
blending also occurs for the P-H 2 O line absorption toward almost all the sight lines surveyed 
here. In these cases, we treated the HF and P-H 2 O emission as features confined to 
the proximity of the background continuum sources and modeled them as an additional 
background continuum emission that adds to the dust emission, as was done in Flagey et ah 
(2013). For these sight lines, we adopted a linear combination of a zeroth-order polynomial 
and up to 3 Gaussian profiles to model the shape of the observed background over the 
regions unaffected by blending with the foreground absorption. In Figs. [T]l6], we plotted our 
best fit models in red over the double sideband spectra of HF and P-H 2 O (black lines) for 
the velocity regions that can be constrained by our data. Since the H 2 O emission profiles 
were constrained using the 1113 GHz features alone, our continuum emission models are 
not as robust as those obtained by Flagey et al. (2013). As a result, the velocity ranges 
significantly affected by these continuum emission features are not discussed further here 
and are not displayed in the bottom two panels of Figs. [Tll6l 

Assuming a sideband gain ratio equal to unity, the flux normalized with respect to the 
continuum flux in a single sideband is then expressed as 


exp(-r) = [TA/T2i(cont) - 0.5]/[Te/TA(cont) 0.5], 
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where Tg is the continuum emission arising from the HF or P-H 2 O -containing gas proximate 
to the star-forming regions and corresponds to the absorption due to the diffuse ISM 
foreground gas of interest for this study. The middle panels in Figs. [T]in] display the resulting 
single sideband, continuum normalized spectrum for HF (black line) and P-H 2 O (blue line) 
versus Vlsr for our sample. The horizontal black lines represent the continuum temperature 
normalized to unity and the zero flux level. One can see that toward all sources, the 
sideband gain ratios are indeed consistent with unity within our uncertainties (see Table |2]). 


3.2. Column Density Measurements 


The HF and P-H 2 O transitions probed here have spontaneous radiative decay rates of 


2.41x10 ^ s ^ and 1.84x10 ^ s respectively. These large rates require high gas densities 
in order for the collisional de-excitation rate to equal the spontaneous radiative decay rate 


for both species. The foreground gas we detect through the HF and P-H 2 O absorptions 
arises from the Galactic diffuse ISM where gas densities have been measured to be at most 
~1000 cm“^ (e.g., Jenkins & Tripp 2001; Sonnentrucker et ah 2007). In the absence of 
a signihcant sub-mm radiation held, we expect that these two species will be entirely in 
their ground rotational states in those foreground gas clouds (e.g., Flagey et ah 2013; 
Emprechtinger et ah 2013). As a result, the optical depth integrated over velocity for a 
gas component detected in absorption via the ground state of HF can be written as (see 


Neufeld et ah 2010) 



Similarly, the optical depth integrated over velocity for a gas component detected in 
absorption via the ground state of P-H 2 O is given by: 


rdn = 4.30 X 10 ^^A^(p — H20/cm^) kms 
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Fig. 1.— Top two panels: Double sideband rms-weighted average spectra for HF and P-H 2 O 
versus Doppler velocity in the Local Standard of Rest frame (Vlsr) toward W3 IRS5 (right) 
and W3(OH) (left). The best £t model to the background continuum plus line emission 
is overlaid in red. The velocity range over which the foreground HF absorption or the 
foreground P-H 2 O absorption is blended with emission proximate to the background source 
is excluded from the modeling. Middle panel: Single sideband continuum normalized spectra 
for HF (black line) and P-H 2 O (blue line). The velocity range in which blending occurs is not 
plotted. Bottom two panels: Optical depth prohles for HF and P-H 2 O (black lines) versus 
Vlsrj respectively. Overlaid in red are the best multi-Gaussian fit models for the HF and 
P-H 2 O prohles for the gas components detected in the foreground of each source that do not 
suffer saturation. 
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Fig. 2.— Same as in Figure 1 for G—0.02—0.07 (+50 km s ^ cloud around SgrA, left) and 
G—0.13—0.08 (+20 km cloud around SgrA, right). 
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G34.3+0.1 W28A 



Fig. 3.— Same as in Figure 1 for G34.3+0.1 (left). In the case of W28A (right), the velocity 
ranges of the foreground absorption and the background emission for HF and P-H 2 O coincide. 
As a result, the optical depth prohles and corresponding multi-component Gaussian hts are 
not presented here. 













































exp(- 


14 


W31C 



p-H,0 l„-0„ ( 1 ^ = 1113.3430 GHz) 


HF 

P-H.0 


-: 

; - n 

Fit 

: / 


L.j 

P-H.0 n 

Fit 11 

^^ll 1 > - 


li 

1 



W33A 



-20 0 20 40 60 

(km s"‘) 


0 20 40 60 80 

(km s"') 


Fig. 4.— Same as in Figure 1 for W31C (left) and W33A (right). 
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G29.96-0.02 W49N 



Fig. 5. 


Same as in Figure 1 for G29.96—0.02 (left) and W49N (right). 
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Fig. 6.— Same as in Figure 1 for W51 (left) and DR21(0H) (right). 
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We used a linear combination of Gaussian components to fit the HF optical depth 
prohles over the LSR velocity ranges where the prohles showed optically thin or moderately 
thick HF absorption depths. For frequencies where the optical depth is large (the absorption 
lines are optically thick), noise fluctuations and small variations in the sideband gain ratio 
result in large uncertainties of the exact values of the optical depth in each frequency bin. 
We place a conservative limit of r = 3 where (2 T^/T^(cont) - 1) < 0.05. For each sight 
line, the gas distribution was modeled by a linear combination of Gaussian components. 
The initial number of Gaussian components was determined by eye and compared to the 
observed spectrum using a y^-minimization technique. After each run, one additional 
Gaussian component was added to the previous model until no signihcant improvement 
to the value was returned by the minimization algorithm. Given the overall similarity 
in component distribution between HF and P-H 2 O, we adopted the best-fit model of the 
HF spectrum (FWHM and velocity range) as our initial guess to model the P-H 2 O optical 
depth prohles for each sight line. We followed the same optimization procedure as for HF 
to derive our best ht model for the H 2 O spectra. We found that for one-half of our sight 
lines, the number of Gaussian components required to best ht the HF and P-H 2 O optical 
depth prohles was identical. The bottom two panels of Figs. [I]l6] display the HF and P-H 2 O 
optical depth prohles (black lines). For clarity, we do not display the portion of each 
spectrum where blending between foreground absorption and background emission occurs. 
Our best ht models are over plotted as red lines only for those absorption components that 
are optically thin or moderately thick. 

While the majority of optically thin or moderately thick absorption features from HF 
and P-H 2 O are htted with a single Gaussian component, some features in our sample require 
the use of a combination of Gaussian components. In the latter case, we sum the column 
densities of the individual components constituting the blended complex and report the 
summed column density over the gas cloud complex. To rehect these diherences in sight line 
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structures, Table [3] reports the LSR velocity range over which the optical depth absorptions 
were integrated rather than the individual Gaussian parameters resulting from modeling 
the sight lines. We only report and discuss measurements for those gas components that 
exhibit optically thin or moderately thick profiles (r < 3) for both HF and H 2 O. 


4. Results 

The spectra displayed in Figs. [T]lH]show that HF is detected in all gas clouds already 
known to trace the diffuse ISM based on HI 21 cm observations (e.g., Lang et ah 2010), or 
on HCO’*' absorption surveys (e.g., Godard et ah 2010). Our survey hence demonstrates 
that the diffuse molecular phase as traced by HF is as ubiquitous as the diffuse atomic phase 
in the Galactic spiral arms. Onr survey also clearly establishes that water is as widespread 
as HF or H 2 in this phase. Toward the Galactic disk sonrces, we hnd that the water velocity 
distribution resulting from onr sight line modeling is identical to that of HF, within 1 km 
s“^, clearly indicating that the water distribntion traces that of HF almost perfectly once 
A^(H20) is in excess of 10^^ cm“^. This striking similarity was noted in earlier work (Nenfeld 
et ah 2010; Sonnentrucker et ah 2010) and is now conhrmed thronghout the Galactic 
volnme we probed. One exception is the gas localized in the immediate vicinity of the 
Galactic Genter where the HF distribntion is observed to be conhned to more discrete gas 
features than the water distribntion (e.g., Lis et ah 2010; Monje et ah 2011, Sonnentrucker 
et ah 2013). For the Galactic disk sources in common between Flagey et ah (2013) and 
this paper, the velocity coincidence of the water absorption featnres with the extent of the 
continnnm emission dne to material close to the sonrce led Flagey et ah (2013) to point 
ont that the water absorption featnres detected toward these sources might not necessarily 
originate in Galactic disk gas. When combining the remarkable similarities in the HF and 
water distributions with the lack of evidence for signihcant HF emission arising from gas 
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local to the background sources, our analysis indicates that the HF and water features we 
considered are mostly tracing the foreground Galactic disk material. Our conclusions are 
consistent with those derived in a recent Galactic [GII] survey by Gerin et ah (2015) who 
showed that the G"*" absorption behaves like a foreground screen to the background sub-mm 
continuum sources. 

Table [3] summarizes our line of sight modeling results and reports on our column density 
measurements as follows. Golumn 1 lists the LSR velocity range over which the optical 
depth prohles for HF and P-H 2 O were integrated. Golumns 2 and 3 report our column 
density measurements for HF, iV(HF), and para-water, iV(p-H20) with corresponding la 
uncertainties in units of 10^^ cm“^. In Golumn 4 we list the total column density of water, 
A^(H20)tot, using an ortho-to-para ratio of 3 for water as measured by Flagey et ah (2013) 
for the Galactic disk. We note that Lis et ah (2010, 2013) measured an ortho-to-para ratio 
around 2.35 for water toward the Galactic Genter source Sgr B2(M). This is the only value 
below the LTE limit that has been found so far and the effect on the H 2 O column density 
is rather minor (~20% decrease). Finally, Golumn 5 reports the abundance of H 2 O relative 
to HF, A^(H20)tot/-^(HF) and associated Icr uncertainty. 

We report column density measurements for both HF and P-H 2 O for a total of 47 
absorption features, all sight lines considered, and limits on either the P-H 2 O or the HF 
column densities for 5 gas components. Thirty of these 47 gas features are satisfactorily 
modeled with a single Gaussian component within our resolution (0.3 km s“^ at the 
HF frequency) and exhibit FWHM ranging from ~1.0 to 4.3 km s“^. The remaining 17 
absorption complexes have ’’effective” FWHM varying from ~4.5 to 10 km s“^. We also 
hnd that the FWHM derived from the P-H 2 O profile modeling are identical to those for 
HF within 1 km s“^, a similarity already noticed in previous work (Neufeld et ah 2010; 
Sonnentrucker et ah 2010) and readily seen in Figs. [T]l6] in this work. Note that this 
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similarity between the HF and H 2 O distributions is genuine throughout the Galactic disk 
and does not result from our line of sight modeling strategy since HF and H 2 O are htted 
separately. Our survey therefore shows that HF and H 2 O trace, in general, the same gas 
volumes throughout the Galactic disk. In the absence of H 2 or GH data, one can thus 
consider using HF to trace H 2 in gas components exhibiting H 2 O absorption, as will be 
discussed in the next sections. 

Figure [7] displays the histogram distribution of the A^(H20)/iV(HF) ratio for the 
47 features that yielded column density measurements for both HF and H 2 O. With this 
sample, we derive a mean water abundance relative to HF of 1.73 with a standard deviation 
around the mean of 0.87. Four velocity ranges exhibit ratios that differ from the mean 
value by more than a factor of 2. When excluding the latter outliers from the distribution 
(gray-shaded histogram), we derive a mean H 2 O to HF abundance ratio of 1.63 with a 
standard deviation around the mean of 0.58. The small differences in the mean values for 
both distributions (less than 12%) indicate that the water abundance distribution relative 
to HF is remarkably narrow throughout the Galactic disk and quite constant to within a 
factor of 2. Gonsidering the non-Gaussian nature of the distribution, we adopted the median 
value of the full (47-point) distribution as our best measurement of the water abundance 
relative to HF, iV(H20)/A^(HF)=1.51, for the Galactic disk. The relative constancy of the 
latter ratio (to within a factor 2.5) was noted earlier by Flagey et al. (2013) for a smaller 
number of gas features. Their measured water abundance relative to H 2 of 5x10“® leads 
to an estimate of the H 2 O/HF ratio of ~1.4, in agreement with our measurements within 
uncertainties. Our survey therefore extends the conclusions of Flagey et al. to a larger 
volume in the Galactic disk. We note that Flagey et al. (2013) used HF/H2=3.6xlO“® 
which will tend to underestimate the H 2 column density and overestimate H 2 O abundance 
relative to H 2 . We also note that the deviations around the median value for the H20-to-HF 
ratio indicate that true variations in the H 2 O/HF ratio do exist. In the following sections. 
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Fig. 7.— Distribution of the H 2 O abundance relative to HF in the Galactic disk gas that we 
sampled. The black-contour histogram displays the distribution using all 47 measurements. 
The mean and median of this distribution are given in the top right corner. Four ratios 
differ from the mean by more than a factor of 2. When excluding those outliers (gray- 
shaded histogram), the distribution of the H 2 O abundance relative to HF yields the mean 
and median values reported in the middle right. 
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we will argue that these variations are mostly due to local variations in the H 2 O/H 2 ratio 
rather than variations in the HF/H 2 ratio in the ISM. As a result, our measurements offer a 
unique opportunity to test the physical processes as well as the chemical pathways involved 
in the production of gas-phase H 2 O in the diffuse interstellar medium. 


5. HF Abundance in the Galactic disk 

In the diffuse ISM, the dehning tracers of the atomic and molecular phases are atomic 
hydrogen (11°) and molecular hydrogen (H 2 ). Direct measurement of II 2 are most easily 
obtained through its ground electronic transitions which he in the far-UV shortward of 
IISOA in the Galaxy. As a resnlt, signihcant effort has been dedicated for decades to 
identifying surrogate tracers of H 2 in the various phases of the ISM. For diffuse molecular 
cloud environments, combinations of FUV and optical observations have demonstrated 
that CH is a valnable tracer of H 2 . The observed CH-H 2 relationship is routinely used to 
estimate H 2 column densities in the absence of II 2 observations. This relationship, however, 
exhibits a large bnt real dispersion of 0.2 dex around its mean value of 3.5x 10-® (Sheffer 
et ah 2008, and references therein; Levrier et al. 2012). The H 2 column densities derived 
from CH measnrements are, thus, only accnrate to within a factor of ~2. The range in H 2 
colnmn densities traced by CH is also limited as the weak CH absorption is typically below 
detection level for diffuse molecular clouds with A^(H 2 ) <fewx 10^° cm“^. CH depletion 
onto dust grains, or “freeze-out” limits the use of CH as a probe of H 2 in denser clouds 
where A^(H 2 ) > 3x 10^^ cm“^ (Mattila 1986). 

Theoretical models of interstellar chemistry have predicted that HF will be the 
dominant reservoir of gas-phase fluorine over a large range of physical conditions (Neufeld 
et al. 2005; Neufeld & Wolhre 2009). In particular, HF is expected to trace H 2 both in 
diffuse gas with molecular fractions of only a few percent (A(H 2 ) >fewx 10^® cm“^, where 
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CH is below the current detection limit) and in dense clouds of larger molecular fractions 
where HF “freeze-out” is not yet signihcant. As a result, HF is predicted to be the most 
sensitive tracer of H 2 in the diffuse molecular regime probed here. One of the goals of 
the PRISMAS key program was to test the latter predictions and determine whether HF 
could be used as a surrogate tracer of H 2 in the sub-mm range where H 2 is not observable 
directly. In the next sections we describe the modifications applied to our diffuse cloud 
models and the updates made to the chemical networks used in the models. We compare 
the HF observations currently available against those new model predictions and we discuss 
our understanding of fluorine chemistry in the diffuse ISM. 


5.1. Model Modifications and Predictions 


We ha ve used a modihed version of the diffuse cloud models presented in iNeufeld et ah 


fennhh and 


grain-surface chemistr y as in 


field attenuation as in 


Neufeld &: Wolhrei (120091). The modifica tions include ice freeze out and 


Hollenbach et ah 


Wolfire et ah 


(120091), the gas-phase reacti ons and radiation 


(120101) , and the oxygen chemistry as in 


Hollenbach et ah 


(120121 ). The model consists of a slab of gas of constant hydrogen nucleus density uh, and 


total width A vtot. wh i ch is i lluminated by the interstellar radiation field, y, measured 


in units of the 


Drainel (119781) held. The gas temperature and the abundances of atomic 


and molecular species are calculated as a function of Ay through the cloud under the 
assumptions of thermal balance, and chemical equilibrium. We assume the slab is embedded 
in an isotropic interstellar radiation held of value x in free-space and thus Y/2_jsjncident 


on op posite sides of the slab. We use the single ray approximation given in 


Wolhre et al.l 


(I 2 OIOI) for attenuation of the isotropic held. 


An important update from the 


Neufeld &: Wolhrei (120091) paper is a new measurement 
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for the rate of the reaction 


F + Ho ^ HF + H 


(ITizniti et ah 


both high or l ow molecu. 
calculations of Zhu et ah 


20141) . This reaction is the dominant formation route for HF in regions with 
a r frac tion. The previously used rate was based on a £t to the 


(120021) for an assumed H 2 population in local thermodynamic 


equilibrium (LTE). The new rate is based on 
and 295 K using a supersonic flow technique. 


experimenta. 


Tizniti et ah 


measurements between 11 K 


(12014!) provide a fit to the rate 


coefficient versus temperature assuming H 2 in LTE. Figure [8] compares the new and old rate 
using the fitted functions. Over temperatures expected in diffuse gas T ~ 50 — 100 K, the 
new rate is a factor of ~ 2.5 times lower than the old rate and thus at low H 2 abundance, 
the new rate results in a lower HF abundance by a factor of ~ 2.5. The HF abundance 
increases with H 2 abundance until all of the gas phase F is locked in HF, at which point the 
HF abundance is insensitive to the H 2 + F reaction rate. 


Additional up da 
and C I with from 


es to the previous codes include fine-structure collision rates for O I 


Abrahamssonetal 


()2nn7l ). and revised photo rates with PAHs using the 
Drainel (|l978l ) radiation field with linear yield functions for the ionization of PAH*^ and electron 


detachment of PAH“. Finally, an important parameter for the gas-phase production of water is 
the cosmic-ray ionization rate, C,. We use the definition that C, is the primary rate per hydrogen 
nucleus. The ionization rate is often quoted as the total rate per H 2 which is a factor ~ 2.3 times 
larger in molecular gas. 


We have run a series of models varying the total extinction through the cloud Ay,tot j the 
density nn, the cosmic-ray ionization rate C,, and the incident radiation field y. The parameters are 


inten 

led to cover a range of values expecte 

in diffuse mo 

ecular clouds (e.g.. 

Sonnentrucker et al 

2007; 

Phillins et al. 

2013; 

GoldsmithI 

2013; 

Indriolo et al. 

2013 

). The cloud extinction varies 


between Ay,tot ~ 0.05 mag (depending on other parameters) and Ay,tot = 4 mag. The density 
varies between nn = 50 cm“^ and nn = 900 cm“^. As a standard model we use a radiation field 
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ig. 8 .— Reaction ra 


:es for HF pr o ducti on using the fitted functions given in iNeufeld et al. 


2 OO 5 I) (old rate) and 


Tizniti et al. 


(12014 f newratel. Bo t 


cloud temperatures in the diffuse ISM, the 


Tizniti et ah 


1 assu me H 2 in LTE. Over typical 


(12014 rate is a factor ~ 2.5 lower 


than the previous rate resulting in lower HF abundances. 


of Y = 1, but examine the case x = 3 as expected for the field in the inner Galaxy (jWolfire et al 
20031 ) . We adopt a standard cosmic-ray ionization rate of ^ = 2 x 10“^® s“^. This rate is close 


to the average rate of ~ 1.5 x 10~^® s“^ estimated from column density measurements over 
50 lines-of-sight by Indriolo &: McCall (2012b). It is also close to the rate of ~ 2.1 x 10“^® 
s“^ estimated by Indriolo et al. (2012a) in the sight line towards W51. Previous diffuse ISM 
studies adopted cosmic ray ionization rates which are lower than the standard value we adopted 
by a factor of ~10. We investigate the effects of lower rates on our predictions as well. For the 


gas-phase abundances of carbon anc 

and 3.9x10"^ ( 

Cartledffe et al. 

2004) 

value of 1.8x 10 

a 

0 

CO 

1 

Neufeld et al. 


oxygen we adopt the values of 1.6x10 ^ (|Sofia et al 


20o4 


20(M). respe ctively. For fluorine, we use a gas-phase abundance 


(j2005l l which amounts to ~ 60% of the solar abundance with 


the remainder depleted onto grains. Our results for the HF abundance will scale linearly with the 
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adopted gas-phase abundance of fluorine. 

Figure 0 shows our model predictions for the HF column density relative to the H 2 column 
density as a function of total visual extinction in the cloud Ay,tot (panel a), as a function of the 
total H 2 column in the cloud iV(H 2 ) (panel b), and as a function of the H 2 fraction in the cloud 
/(H 2 ) = 2 iV(H 2 )/[iV(H°) -|- 2A^(H2)] (panel c). For the models shown in FigO y is hxed at x =1 
in units of the Draine (1978) field. The A^(HF)/A^(H 2 ) ratio is predicted to rise monotonically 
with increasing Ay,tot, ^(^^2), and /(H 2 ) and varies by at most a factor of ~ 3 while C and nn 
vary by at least a factor of 10 over the entire parameter range considered here. 


The dependence of the A^(HF)/A^(H 2 ) ratio on total cloud extinction Ay^tot can be understood 
by considering the dominant formatioi 
detailed HF chemistry is presented in 


Neufeld et al. 

0 

0 

CM 

) and 

Neufeld & Wolfire 

(20091 


the production of HF proceeds by the reaction of H 2 with F, while the destruction occurs through 
reactions with C"*", Hs"*", He"*", and Si"*" plus photodissociation. For the low density models, C'*' 
dominates the destruction of HF, to Ay^t^t ~ 1.5 mag, at which point electrons recombining with 
C"*" lead to a lower 0"*“ abundance and higher HF abundance. At high density (nn > 300 cm“^), 
the recombination of C"*" is drawn to the cloud surface leading to higher HF columns even for 
clouds of low Ay,tot- Increasing cloud column densities results in a drop in the C"'' abundance and 
reactions with less abundant ions dominate the destruction of HF. Thus the HF destruction rates 
fall and the A^(HF)/A'(H 2 ) ratio rises. 


Since both the HF column density and H 2 column density depend on the H 2 abundance, the 
ratio does not depend directly on the H 2 abundance. In addition, since both the formation of HF 
and destruction of HF have the same dependence on density (oc n^), the HF fractional abundance 
is not directly density dependent. As a result, we only see a weak density dependence in Figure [H 


The HF abundance predictions are also mostly insensitive to variations in Note that a 
slight increase of the HF abundance of at most 10% is predicted at Av',tot= 4 mag for ran < 100 
cm“^ when increasing by a factor of 10. Cosmic rays play a role in the destruction of both HF 
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Fig. 9.— Predictions for the HF abundance relative to H 2 in a single cloud as a function 
of total visual extinction {Ay,tot, panel a), total H 2 column density (A^(H 2 ), panel b), and 
molecular fraction (/(H 2 ), panel c) using our modified 2-sided PDR model. The UV held 
is given in units of the Draine (1978) held, is hxed to y =1 and is assumed to be incident 
isotropically at a value of y = 1/2 on each side of the cloud. In panel a), we list the range 
in cosmic-ray ionization rate ( (10“^® s“^) and density tt-h (cm“^) we explored. In panel c), 
colored symbols show model predictions for clouds with Ay, tot = 0.1, 0.2, 0.3, 0.5, 1.0, 1.5, 
2.0, 2.8 and 4.0 mag. For models (uh = 300 cm“^ and C, = 0.7x10“^® s“^) or (uh = 900 
cm“^ and C = 2x10“^® Ay^t^t = 0.05 mag is also shown. To convert Rytot to hydrogen 
nuclei column density iVn = A^(H°) -|- 2 A^(H 2 ) we use iVH(cm“^)= 2x10^^ ^ytot- 









and H 2 , as they produce H 3 + and He"*" which then react with HF and H 2 . The decrease in the H 2 
abundance is greater than that in the HF abundance which, in turn, results in the ~10% increase 
in the iV(HF)/A^(H 2 ) ratio seen in Fig. [9l 


5.2. HF Observations 

The grid of model calculations presented in the previous section predicts that the column 
density of HF with respect to H 2 varies between A^(HF)/A^(H 2 ) ~ 0.9 xl0“® at the cloud surface 
{Ay,tot ~ 0.08 mag) in the low density regime to A^(HF)/A^(H 2 ) ~ 3.3 x 10“® at the cloud center 
{Ay,tot ~ 4 mag) in the high density regime considered here. As a result, HF is predicted to be a 
very sensitive probe of H 2 , as it can trace H 2 both in mostly diffuse atomic gas (/(H 2 )< 0 . 1 ) and 
in mostly diffuse molecular gas (0.1< /(H 2 ) <0.9) in the interstellar medium. The HF abundance 
in our models never reaches the value of 3.6 x 10“® expected if all gas-phase F is locked into 
gas-phase HF; a value derived from FUV measurements of the abundance of gas-phase atomic 
fluorine (e.g.. Snow et al. 2007). The leveling off around a value of 3.3 xl0“® seen in Fig. [9]for nn 
= 900 cm“^, is due to the incomplete conversion of F into HF in the outer regions of the cloud. If 
the gas becomes fully molecular slightly before all F is locked in HF, then the total column density 
ratio A^(HF)/A^(H 2 ) is less than the F abundance. 

Following the first Herschel/mFl detections of the ground rotational transition of HF 
toward sight lines comprised mostly of diffuse atomic/molecular clouds (Neufeld et al. 2010; 
Sonnentrucker et al. 2010), numerous additional HF detections were reported both with the high 
resolution HIFI data and the lower resolution SPIRE and PACS data (e.g., Kirk et al. 2010). 

At this point in time, HF has been detected in the Galactic Center (e.g., Sonnentrucker et al. 
2013; Goicoechea et al. 2013), in hot cores around massive proto-stars (e.g., Phillips et al. 2010; 
Emprechtinger et al. 2012), in the ISM of extra-galactic sources up to z=2.6 (e.g., Monje et al. 
2011; Monje et al. 2014) and in the Galactic disk (Flagey et al. 2013; this work). These studies 
clearly established the ubiquitous presence of the HE molecule in the interstellar medium, as 
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predicted by our models and as expected from a tracer of H 2 . 

The A^(HF)/A^(H 2 ) values currently available in the literature range from 1.0 x 10“® to 2.5x 
10“® (e.g. Sonnentrucker et ah 2010; Phillips et al. 2010, Monje et ah 2011; Emprechtinger 
et al. 2012) and are fully consistent with our model predictions. With the exception of the 
HE/H 2 ratio estimated by Indriolo et al. (2013) and discussed below, those ratios were estimated 
indirectly from simultaneous measurements of HE and either CH or along these sight lines. 

The CH column densities were converted into H 2 columns using the most recent estimate of the 
CH-H 2 relationship A^(CH)/A^(H2)=3.6 x 10“® (Sheffer et al. 2008). As mentioned earlier, the 
CH-H 2 relationship exhibits a standard deviation of about 0.2 dex that is real and thought to 
be related to the chemical pathways involved in the formation of CH in the diffuse ISM. The 
^^CO column densities where converted to H 2 columns for a given carbon Galactic isotopic ratio 
and for a given CO-to-H 2 ratio; the latter two quantities also vary by factors of a few in the 
Galactic disk (Sonnentrucker et al. 2007; hanger & Penzias 1990). With these caveats in mind, it 
is interesting to note that the range in the A^(HF)/Ai(H 2 ) ratio measured to-date with Herschel 
is fully consistent with our model calculations, thus, validating our understanding of both the 
physical conditions and the new reaction rates involved in the chemistry of interstellar fluorine. 

A few direct measurements of A^(HF)/A^(H 2 ) were recently obtained by Indriolo et al. (2013) 
for a handful of targets sampling diffuse molecular clouds where A^(H 2 ) and A^(HF) were derived 
from FUV and near-IR observations, respectively. For all sight lines where measurements of both 
quantities were obtained, Indriolo et al. (2013) derived A^(HF)/A^(H 2 ) ratios between 0.5 and 1.4 
X 10“®. Of particular interest are the results obtained toward the translucent sight line HD154368. 
This sight line was best modeled to be comprised of diffuse molecular material with nH= 325 
cm“^, x=3 and (= 0.5x 10“^® s“^ by Spaans et al. (1998) while FUV measurements based on C 2 
observations yielded nH= 240 cm“^ (Sonnentrucker et al. 2007). Our model calculations for the 
HF abundance using the parameters from Spaans et al. (1998) and an average density of nn = 300 
cm“^ yield a ratio of A^(HF)/Ai(H 2 ) = 2.14 x 10“® toward this sight line. Our model predictions 
are consistent with the direct measurements of Indriolo et al. (2013) of A^(HF)/V(H 2 ) = 1.15 
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± 0.41 within 3 -cj. In summary, all measurements obtained to date -whether direct or indirect- 
corroborate the prediction that HF can be used as a surrogate tracer of H 2 in the diffuse ISM, 
down to very low molecular fractions where CH is typically below our detection level. Since direct 
H 2 measurements are not available for the entire set of HF and H 2 O absorption features we report 
in this survey, we will use the HF column densities we derive as proxies for H 2 in the remainder of 
the paper. 


6. H 2 O Abundance in the Galactic disk 


6.1. Model Predictions 


Figure do] shows the predictions of the PDR model described in the previous section for the 
H 2 O column density relative to the H 2 column density as a function of total visual extinction in 
the cloud {Ay,tot, panel a), as a function of the total H 2 column density in the cloud (A(H 2 ), 
panel b) and as a function of the molecular fraction in the cloud (/(H 2 ), panel c). The UV field is 
given in units of the Draine (1978) field and is fixed to x=l. 


T he Ho O che mistry was discusser 


(120091 . 


2013 ), and 


van Dishoeck et al 


in det ail by 


van Dishoeck &: Black (1986), 


Hollenbach et al 


(|2013l b As in the case of HF, the variations in the 


H 2 O abundance reflect the competition between the formation and destruction processes that 
contribute over the parameter range explored here. In brief, the destruction of H 2 O is dominated 
by photodissociation and, at low Ay, has some contribution by reaction with C+. The production 
of H 2 O proceeds by a combination of surface chemistry on grains and ion-neutral chemistry in 
the gas phase. When the molecular fraction is low the ion-neutral chemistry in the gas phase 
is initiated by cosmic-ray ionization of atomic H, while when the molecular fraction is high, 
the ion-neutral chemistry is initiated by cosmic-ray ionization of H 2 . In both the atomic and 
molecular branches the rate of production of H 2 O increases with the H 2 abundance since both 
branches require reactions with H 2 to proceed in the gas phase. When both the H 2 and electron 
abundances are high, the production of water can be reduced due to electron recombinations 
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Fig. 10.— Predictions for the H 2 O abundance relative to the H 2 in a single cloud as a 
function of total visual extinction (Aytot, panel a), total H 2 column density (A^(H 2 ), panel 
b), and molecular fraction (/(H 2 ), panel c) using our modihed 2-sided PDR model. The 
UV held is given in units of the Draine (1978) held, is hxed to y =1 and is assumed to be 
incident isotropically at a value of x = 1/2 on each side of the cloud. In panel a), we list the 
range in cosmic-ray ionization rate ( (10“^® s“^) and density tt-h (cm“^) that we explored. 
In panel c), colored symbols show model predictions for clouds with Ay, tot = 0.1, 0.2, 0.3, 
0.5, 1.0, 1.5, 2.0, 2.8 and 4.0 mag. For models (nn = 300 cm“^ and C = 0.7x10“^® s“^) or 
(riH = 900 cm“^ and C = 2x10“^® Ay^t^t = 0.05 mag is also shown. To convert Rytot to 
hydrogen nuclei column density Nu = iV(H®) -|- 2At(H2) we use iVH(cm“^)= 2x10^^ Rytot- 
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with Hg . When the H 2 abundance is low, but the electron abundance remains high, the water 
production can be reduced by electron recombination with H'*' or by neutralization of by 
PAH“ or by PAH. Finally, when the H 2 abundance is high, but the electron abundance is low, the 
gas-phase production of H 2 O proceeds efficiently through ion-molecule reactions. 

To further explain the physical regimes where the various processes dominate, we start with a 
typical diffuse cloud visual extinction (Ay^tot = 0.5 mag), density (nn =100 cm“^), and cosmic-ray 
ionization rate (C = 2x10“^® and then discuss the variation in the H 2 O production, while 

varying these parameters. For the Ay^tot ~ 0.5 mag and nn < 100 cm“^ models, the production of 
H 2 O is dominated by ion-neutral chemistry through the cosmic-ray ionization of atomic hydrogen. 
To first order the H 2 O abundance scales as C/^H and thus increasing density results in a lower 
A^(H 20 )/A’(H 2 ) ratio, for a given cosmic-ray rate. 

For models u-h > 300 cm“^, the molecular fraction rises and the atomic hydrogen abundance 
falls sufficiently so that the ion-molecule chemistry proceeds mainly through cosmic-ray ionization 
of H 2 . However, for these models, the electron abundance remains high and the gas-phase 
production of H 2 O is reduced due to the recombination of Hjj", an intermediary in this particular 
H 2 O production route. As a result, H 2 O is produced mainly by reactions on grains. In the limit of 
H 2 O produced by grain chemistry alone and destroyed by photodissociation, the A'(H 20 )/A'(H 2 ) 
ratio is a function of nn/x (Hollenbach et ah, 2009). 

For the Ay^tot ~0.5 mag, nn =300 cm“^, and C =0.7x10“^® s“^ model, only ~ 15% of the 
water is predicted to be produced through gas-phase reactions. For the lowest cosmic ray rate 
we tested (C = 0.2 x 10“^® s“^), the models with u-h > 100 cm“^ and x = 1 dominated by 
grain surface chemistry and increasing the density nn produces more water and, hence, higher 
iV(H20)/iV(H2). 

At lower Ay^tot (< 0.5 mag), curves for which grain surface chemistry dominates are seen to 
rise. Since the grain surface chemistry reactions do not depend on the H 2 abundance, the curves 
rise as the H 2 abundance drops. At higher Ay^tot > 2 mag, for all models, the electron abundance 
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drops sufficiently due to absorption of FUV photons, so that water can be produced efficiently by 
ion-molecule chemistry. 

We note that Dulieu et al. (2013) have presented experimental evidence to suggest that 
OH and H 2 O formed on bare silicate grain surfaces can be chemically desorbed upon formation. 
We consider the effects of chemical desorption by assuming that 30% of the OH is immediately 
desorbed and 70% of the H 2 O is immediately desorbed (see their Fig. 4). We find at most a ~ 
60% increase in the H 2 O column density over our model parameter space. The effect is largest 
where surface chemistry dominates the production of H 2 O and for the smallest cloud column 
densities. The increase drops to ~ 40% for a total cloud column density of Ay,tot =2.0, and drops 
to ~ 10% for a total cloud column density of Ay,tot = 4.0. We also note that Minissale et al. 
(2013), and Minissale et al. (2014) have found that 0-atom diffusion can be quite rapid on grain 
surfaces leading to the production of O 2 and O 3 . This process is important deep in molecular 
clouds where the incidence of oxygen atoms exceeds that of H. These conditions are, however, not 
applicable to the models presented here. 


6.2. H 2 O Observations 

Since H 2 is not observable directly in the sub-mm, efforts are made to find surrogate tracers 
of H 2 in this wavelength domain in order to obtain abundance measurements relevant to test our 
knowledge of diffuse cloud chemistry. In Section 5, we demonstrated that HF can be used as a 
surrogate tracer of H 2 , down to much lower molecular fractions than the more typical tracers such 
as CH. Since we do not have H 2 or CH column density measurements for all 47 features presented 
here, we use HF as a tracer of H 2 and discuss the H 2 O abundance relative to HF -A^(H 20 )/A^(HF)- 
in the remainder of this paper. 

As described in Section 4, we detected absorption from H 2 O toward all 12 sight lines included 
in our survey and we measured its column density in a total of 47 cloud components (see Table [3]) . 
We compared our H 2 O column density measurements against those reported by Flagey et al. 
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(2013) for the 17 gas components in common to both studies; these are detected in the foreground 
to W49N, W51, G34.3+0.1, W33A and DR21(OH). We find that the two sets of measurements 
are consistent within 2 a for 13 components. For four components, the H 2 O column densities we 
measured are either lower by 50% (DR21(OH)) or higher by up to a factor of 3 (G34.3+0.1 and 
W49N) than those reported by Flagey et al. (2013). In all four cases, the background continua 
local to these gas components are somewhat affected by the presence of water emission local to 
the sources; the discrepancies we noted above are readily explained by the different methods we 
adopted to account for the background emission toward each source and give a measure of the 
modeling uncertainty in such cases. The H 2 O column density measurements we report for these 4 
velocity components (marked with an asterisk in Table [3|) should be used with caution. 

Flagey et al. (2013) reported no dependence of the H 2 O abundance relative to H 2 on 
Galactocentric distance for the 6 sources they considered. For those sources located in the inner 
Galaxy in our extended survey, we computed a Galactocentric distance for each cloud we detected 
in HF and H 2 O based on the source longitude and Vlsr to see whether the variations in the 
A^(H20)/At(HF) we measure show a dependence on position in the Galaxy. We find no dependence 
of the A^(H20)/At(HF) ratio with Galactic radius either, a conclusion consistent with that of the 
smaller Flagey et al. (2013) survey. Our results therefore add weight to the conclusion that the 
variations in the A^(H 20 )/A^(HF) ratio we measure can be considered representative throughout 
the Galactic disk. With a median value of A^(H20)/A^(HF) = 1.51 in diffuse clouds, H 2 O also 
appears as an alternative tracer of H 2 (within a factor of 2.5) in the absence of HF or CH spectra, 
as already suggested in Flagey et al. (2013). 

The comparison of our model predictions displayed in Fig. [9] and Fig. [TOl shows that, in the 
diffuse molecular cloud regime (/hj < 0.5), the column density ratio of H20-to-HF varies from 
0.5 < A^(H20)/A^(HF) < 3, all models considered. Our models predict that the water abundance 
relative to HF varies by up to a factor 6, depending on the local gas physical conditions. The 
model comparisons further show that the variations in the A^(H 20 )/A^(H 2 ) ratio are much larger 
(factor of 2 to 4) than the variations in the A^(HF)/A^(H 2 ) ratio (at most factor of 2) both as a 


35 


function of cloud column density for a given set of parameters and between models with different 
inpnt parameters. Our models further predict that the variations in the A^(H20)/A^(HF) are 
mostly driven by variations in the H 2 O column density, rather than in the HF column densities 
for a given set of model parameters. 

In Figure fm we display the iV(H20)/A^(HF) ratio versus iV(HF) for all 47 foreground clouds 
with optically thin or moderately thick absorptions. The measurements toward each sight line 
are represented by a particular symbol and color that are listed at the top of the panel. We 
over-plot our model predictions for the A^(H 20 )/iV(HF) ratio onto our measurements using the 
same symbol and color coding as in FigsiQl&llOl The range of model parameters we use is listed at 
the top left of the panel. We use a log-log scale for both axes for clarity purposes, la uncertainties 
are shown for all measurements and 3a uncertainties are plotted as downward arrows for upper 
limits in the water abundance. This figure shows that the iV(H20)/A^(HF) ratio varies by about a 
factor 5 overall; our observational results are therefore fully consistent with the model predictions 
presented here and known to best bracket the physical conditions in diffuse clouds. 

Our sight line analyses (Section 2) indicated that HF and H 2 O are co-located in the gas 
components we detect throughout the Galactic disk, meaning that both molecnles are subject 
to the same physical conditions within each gas clump along the sight lines. As a result, the 
variations we measure are not caused by differences in the spatial distributions of HF or H 2 O 
along a given sight line (geometric effects) but are caused by true variations in the molecular 
content for a given gas component. 

In particular, Fig. [11] displays a subtle break in the water abundance relative to HF 
once A^(HF)=5xlO^^ cm~^. For iV(HF) < 5x10^^ cm“^, our measurements indicate that the 
A^(H20)/A^(HF) ratio varies by about a factor of 6 to 9 from ~ 0.5 to 4.6. The lowest water 
abundances relative to HF are best reproduced by models with low cosmic ray ionization rates 
(C ^ 0.7x10“^® s“^) over the density range known to be representative of the diffuse ISM while 
the highest water abundances relative to HF are best represented by ~ 2x10 s ^ and uh < 
100 cm“^. The prevalence of mixtures of physical conditions within the galactic disk diffuse ISM 
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Fig. 11.— H 2 O column density relative to that of HF versus HF column density throughout 
the Galactic disk probed by our Herschel survey. Each sight line is represented by a particular 
symbol and color which are listed at the top of this hgure. Overplotted on our measurements 
are the predictions obtained from a series of 2-sided PDR model runs. The parameter space 
explored in our runs is listed in the upper left of the panel with ( in units of 10“^® s“^ and uh 
in units of cm“^. 1 a uncertainties are over plotted for each measurement. S-a uncertainties 
are displayed for all upper limits. The colored vertical bars over plotted on each model run 
correspond to a depth in the cloud of Ay,tot ~ 0.5 mag. The upper limit for HD154368 was 
derived using H 2 O results from Spaans et ah (1998) and HF and H 2 results from Indriolo et 
ah (2013). 
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and along any given sightline has long been recognized and is clearly evidenced here again, as 
multiple model runs can accommodate our measurements for this particular HF column density 
range. 

For A^(HF) > 5x10^^ cm“^, however, our measurements are not consistent with models 
using ( < 0.7x10“^® s“^. As a matter of fact, over 80% of our measurements in this specific HF 
column density range show water abundances relative to HF consistent with model predictions for 
single clouds with ( = 2x10“^® s“^, y =1 and uh <100 cm“^. The variations in the measured 
A^(H20)/A^(HF) in this HF regime amount to at most a factor of 3, thus, pointing toward a 
population of low gas density clouds with properties intrinsically very similar. Considering that 
the HF abundance relative to H 2 currently measured varies by at most a factor of ~2 (see Section 
5), we conclude that the variations we measured in the A^(H20)/A^(HF) ratio are mostly driven 
by variations in the H 2 O/H 2 ratio within the gas clouds we probe. 

7. Discussion 

In the diffuse ISM, the production of H 2 O within our models proceeds via both gas-phase 
reactions and grain-surface chemistry. Since the contribution of each mechanism depends on the 
cloud physical conditions, our measurements give us a unique opportunity to test the chemical 
pathways predicted to lead to H 2 O formation in this regime. 

In Figure dH we show the fractional rate of H 2 O production on grain surfaces compared to 
the total production rate of H 2 O as a function of A^(HF). The rates are evaluated at cloud center 
where the H 2 O abundance is generally highest and dominates the contribution to the column 
density. These curves can be directly compared with those in Figure [TTJ 

As already mentioned in Section 6, at a column of A^(HF) ~ 5x10^^ cm“^ and for nn < 
100 cm“^, the production of H 2 O is dominated by ion-neutral chemistry (see Fig |I2p . The H 2 O 
abundance scales as C/^H and thus a larger density results in a lower A^(H 20 )/A^(HF) ratio for 
a given cosmic-ray rate (Fig. [TT]l . At higher densities nn > 300 cm“^, grain-surface reactions 


38 


dominate the production of water. Figure [12] shows that for riu = 300 cm ^ and C = 0.7x10 
s“^ only about 15% of the water is produced through gas-phase reactions. 

For A^(HF) < 5 x 10^^ cm“^ (^v,tot < 0.5 mag) , the curves for which grain surface chemistry 
dominates are seen to rise in Fig. 10. Since the grain chemistry reactions do not depend on the H 2 
abundance while the HF production does, the curves rise as the H 2 and thus HF abundance drops. 
At higher Ay,tot) the electron abundance drops sufficiently due to absorption of FUV photons, so 
that water can then be produced efficiently through the cosmic-ray ionization of H 2 . The curves 
thus rise due to an increasing abundance of H 2 O. 

Note that in the limit in which H 2 O is produced by grain surface chemistry alone and 
destroyed by photodissociation, while HF is both formed and destroyed by gas-phase reactions, 
the ratio A^(H20)/A^(HF) is a function of nn/x (see Figs. [IT] fc: [T^ . For the lowest cosmic ray 
rate we considered here (C = 0.2 x 10“^® s“^), the curve for tt-h = 100 cm“^ and x = 1 lies close to 
the curve for ^ = 0.7 x 10“^® s“^, nn = 300 cm“^ and x = 3, as expected for H 2 O production by 
grain chemistry (see Fig. m- 

Our measured water abundances relative to HF reported in Table [3] are all bracketed by model 
predictions based on a grid of gas physical conditions known to be representative of the diffuse 
ISM. At low visual extinctions and for densities greater than tt-h =100 cm“^, our models predict 
that grain-surface chemistry is necessary to account for some of our observations. For Ay,tot > 0.5 
mag, our models indicate that contributions from both gas-phase and grain-surface chemistry are 
required to account for our measurements as well. The contribution of each production mechanism 
depends intimately on the electron density distribution within the cloud, as recombination with 
H+ and H^^ inhibit the gas-phase production routes. 

Within our uncertainties, most measurements are compatible with more than one model. 
These degeneracies are mainly due to two factors. Firstly, some of our model predictions are 
degenerate and measurements of HF and H 2 O alone do not allow us to lift the model degeneracies; 
additional observational probes such as HCO+ (e.g., Godard et al. 2014) are necessary to do so. 
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Fig. 12.— Fractional rate of H 2 O production on grain surfaces compared to the total pro¬ 
duction rate of H 2 O versus HF column density in the column density range probed by the 
survey. The parameter space explored in our runs is listed in the upper section of the panel 
with ( in units of 10“^® s“^ and uh in units of cm“^. The UV held is given in units of the 
Draine (1978) held and is hxed to y =1 or y =3. 
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Secondly, our measurements are not obtained toward true single clouds contrary to the model 
assumptions, as discussed earlier. To compare the water abundance relative to HF versus HF 
column density on equal footing one would have to determine the HF column density in each 
gas clump along a given sight line, a task which is observationally very challenging. One rare 
exception might be the sight line toward HD154368. As mentioned earlier, Spaans et al. (1998) 
and Rachford et al. (2002) determined that true translucent conditions are present in the gas 
intercepting the sight line to this star. The decomposition performed by Rachford et al. (2002) 
suggests that ~73% of the H 2 gas observed with FUSE data belongs to the translucent phase. 

As a result, it is reasonable to assume that ~73% of the total HF column density we measure is 
also associated with this translucent phase. Correcting the total measured HF column density 
by that fraction has the effect of moving the upper limit on the A^(H 20 )/A^(HF) ratio to the 
left (or A^(HF)~ 1.2 xlO^^ cm“^) in the bottom panel of Fig. fTTl One can see that the overall 
impact of such a sight line “decomposition” on our measurements would primarily be to shift the 
distribution horizontally and to the left - or to lower A^(HF). 

Despite these degeneracies, it is important to note that the lower envelope on the 
A^(H 20 )/A^(HF) ratios we measure is too high to understand without grain-surface chemistry 
being important. The overall range in A"(H20)/A^(HF) ratios that we observe in itself demonstrates 
the importance of grain surface chemistry in the production of water in this diffuse ISM. It is also 
remarkable to see that the adopted physical conditions and chemical network that best bracket our 
observations remain largely unaltered, when taking the latter distribution shift into account. Our 
measurements therefore indicate that the range in densities, FUV and cosmic-ray ionization rates 
required to model our observations are not only pervasive throughout the diffuse gas composing 
the Galactic disk but they are also quite tightly constrained. Finally, our results also validate our 
understanding of the various chemical pathways thought to lead to the gas-phase abundance of 
HF and H 2 O in the diffuse ISM. While most chemical pathways involve gas phase reactions in 
diffuse clouds, our measurements do confirm that grain surface chemistry also plays a significant 
role in the production of gas-phase water in such low density environments. 
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We note that the effects of turbulent dissipation regions (TDR) on the chemistry in diffuse 
clouds were recently extensively tested and were found to successfully reproduce the wealth of CH''” 
and SH"*" absorption measurements obtained with Herschel toward sight lines probing the diffuse 
ISM (Godard et al. 2014 and references therein). The production of gas-phase water was also 
included in those models. Godard et al. (2014) found that for gas densities in the range we probe 
here (nH= 50-300 cm“^) and for ^ytot~0.5 mag and a primary cosmic ray rate per H of 0.4x10“^® 
s“^ cm“^, TDR could contribute up to 50% of the total (TDR plus ion-molecule) gas-phase 
production of water. The Godard et al. (2014) models do not include the effect of grain-surface 
chemistry, while our models do not include the effects of TDR on the water production, hence a 
direct comparison of these model predictions is not feasible. However, all differences considered 
and for this particular set of parameters, our models predict that grain surface chemistry produces 
up to 85% of the observed gas-phase water abundance, with ion-molecule reactions producing 15% 
of the observed gas-phase water abundance. Consequently, ion-molecule reactions combined with 
TDR effects alone (without any contribution from grain surface chemistry) would fall short (~35% 
of that required) in terms of reproducing the observed water abundance measured in this survey. 


8. Conclusions 

We have presented the most comprehensive survey of HF and H 2 O observations obtained 
in the diffuse interstellar medium of the Galactic disk using the Herschel/HIFl instrument. The 
column density of both molecules is measured toward 47 discrete gas components detected in 
absorption and probing the Galactic disk volume. We demonstrate that both molecules are 
ubiquitous in the diffuse gas of our Galaxy. We find that the HF and H 2 O velocity distributions 
trace each other almost perfectly, in the disk, establishing that HF and H 2 O essentially probe the 
same gas-phase volume. 

We compare our observations to state-of-the-art diffuse cloud models that were modified to 
include grain surface chemistry, as well as updated reaction rates for the fluorine chemistry. Our 
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measurements corroborate theoretical predictions that HF is a very sensitive tracer of H 2 down to 
molecular fractions of a few percent, a regime where more common tracers such as CH are below 
detection level. 

We have used HF as a surrogate tracer of H 2 to study the variation of the H 2 O column 
density -relative to HF- within the Galactic disk diffuse gas. We find that the A^(H20)/iV(HF) 
ratio shows a narrow distribution with a median water abundance relative to HF of 1.51. Our 
results therefore add weight to the previous suggestion that H 2 O can also be used as tracer of H 2 - 
within a factor of 2.5- in the diffuse ISM, in the absence of HF or CH observations. 

We have demonstrated that the overall variations of a factor of 2.5 around the median in the 
A^(H 20 )/A^(HF) ratio are driven by true variations in the H 2 O column density within the disk. 
Our measurements, therefore, provide us with a unique opportunity to test the chemical pathways 
predicted to play a role in the water production in the diffuse ISM. We find that the range in water 
abundances relative to HF that we measure could only be explained if significant grain surface 
chemistry production occurred in addition to gas-phase ion-molecule production. While most 
chemical pathways involve gas phase reactions alone in the diffuse ISM, our survey confirms that 
grain surface chemistry can play a significant role in the production of some molecular species, 
such as gas phase H 2 O, in this low density environment. 
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Table 1: Target List. 


Target 

(.) 

RA 

n (') 

(°) 

Dec 

nr) 

1 

n 

b 

n 

D 

(kpc) 

Ref. 

W3 IRS5 

02 

25 

40.6 

+62 

05 

51.0 

133.715 

+1.215 

1.83 

1 

W3(OH) 

02 

27 

03.8 

+61 

52 

25.0 

133.946 

+1.06 

2.04 

2 

G-0.02-0.07 

17 

45 

50.2 

-28 

59 

53.0 

359.97 

-0.07 

8.34 

3 

G-0.13-0.08 

17 

45 

37.4 

-29 

05 

40.0 

359.30 

-0.432 

8.34 

3 

G34.3+0.1 

18 

53 

18.7 

+01 

14 

58.0 

34.26 

+0.15 

3.80 

4 

W28A 

18 

00 

30.4 

-24 

04 

00.0 

5.9 

-0.39 

1.28 

5 

W31G 

18 

10 

28.7 

-19 

55 

50.0 

10.62 

-0.38 

4.95 

6 

W33A 

18 

14 

39.4 

-17 

52 

00.0 

12.91 

-0.26 

2.40 

7 

G29.96-0.02 

18 

46 

03.9 

-02 

39 

21.9 

29.96 

-0.02 

5.26 

8 

W49N 

19 

10 

13.2 

+09 

06 

12.0 

43.17 

+0.01 

11.1 

9 

W51 

19 

23 

43.9 

+14 

30 

31.0 

49.49 

-0.39 

5.41 

10 

DR21(OH) 

20 

39 

01.0 

+42 

22 

48.0 

81.72 

+0.57 

1.50 

11 


References: 1- Imai et al. (2000); 2- Hachisuka et al. (2006) ; 3- Reid et al. (2014); 4- Fish 
et al. (2003); 5- Motogi et al. (2011); 6- Sanna et al. (2014); 7- Immer et al. (2013); 8- 
Zhang et al. (2014); 9- Zhang et al. (2013) ; 10- Sato et al. (2010); 11- Rygl et al. (2012) 








Table 2. Observation Summary. 


Target 

tea:p(HF) 

(s) 

TA(cont)[HF]“ 

(K) 

texp(H20) 

(S) 

TA(cont)[H2 0]“ 

(K) 

Observation Dates 

HF & H 2 O Obs IDs 

starting with 1342: 



W3 IRS5 

492 

3.56 

± 

0.10 

2965 

3.11 

± 

0.02 

2012-07-20 

& 

2010-07-18 

248382-83-84 & 201591 

W3(OH) 

9275 

4.46 

± 

0.03 

3540 

4.14 

± 

0.03 

2013-03-27 

& 

28 

268474-75-76 & 268605-06-07 

G-0.02-0.07 

3693 

0.93 

± 

0.04 

2139 

0.92 

± 

0.04 

2010-10-06 



205885-86-87 & 205882-83-84 

G-0.13-0.08 

746 

1.43 

± 

0.11 

428 

1.39 

± 

0.11 

2011-09-14 



228613-14 & 228615-16 

G34.3+0.1 

522 

9.20 

± 

0.12 

246 

7.57 

± 

0.10 

2010-04-18 



195074-75-76 & 195070-71-72 

W28A 

562 

6.03 

± 

0.10 

426 

4.71 

± 

0.06 

2011-03-11 



215860-61-62 & 215863-64-65 

W31G 

534 

7.74 

± 

0.11 

306 

5.99 

± 

0.09 

2010-03-05 



191690-91-92 & 191687-88-89 

W33A 

3775 

2.36 

± 

0.04 

2694 

1.91 

± 

0.03 

2011-03-11 



215875-76-77 & 215872-73-74 

G29.96-0.02 

9162 

3.20 

± 

0.03 

2824 

2.69 

± 

0.01 

2013-03-27 

& 

2011-09-29 

268483-84-85 & 229875-76 

W49N 

522 

11.48 

± 0.18 

246 

9.06 

± 

0.16 

2010-03-22 



192595-96-97 & 192592-93-94 

W51 

9259 

11.12 ± 0.04 

3582 

10.51 

± 0.04 

2013-03-27 

& 

28 

268478-79-80 & 268611-12-13 

DR21(OH) 

1014 

5.36 ± 0.06 

1086 

4.56 ± 0.05 

2010-05-12 



196502-03-04 &196506-07-08 


T^(cont) is the double sideband continuum antenna temperature. 
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Table 3. Column Density Measurements. 


Dlsr Range 

(km s~^) 

iV(HF) 

(10l2 cm-2) 

iV(p-H20) 

(10^^ cm“^) 

iV(H20)tot 

(10^^ cm“2) 

Af(H20)tot 

Ar(HF) 

W49N 

+30,+34* 

12.74+ 0.22 

4.82 + 0.16 

19.27 + 0.66 

1.51 + 0.06 

+43,+48 

2.91 ± 0.25 

1.47 + 0.27 

5.87 + 1.07 

2.02 + 0.41 

+49,+53* 

5.73 ± 0.27 

1.80 + 0.18 

7.19 + 0.73 

1.25 + 0.14 

+53,+55 

6.37 ± 0.24 

3.04 + 0.25 

12.18 + 1.01 

1.91 + 0.17 

+55,+58 

5.63 ± 0.22 

3.80 + 0.28 

15.22 + 1.11 

2.70 + 0.23 

+64,+71 

8.36 ± 0.89 

2.82 + 0.88 

11.27 + 3.51 

1.35 + 0.44 

W51 

+4,+6 

4.88 ± 0.17 

0.88 + 0.16 

3.53 + 0.62 

0.72 + 0.13 

+6,+10 

11.39+ 1.07 

5.37 + 0.11 

21.49 + 4.54 

1.89 + 0.44 

+11,+15 

1.46 + 0.83 

0.22 + 0.08 

0.88 + 0.32 

0.60 + 0.40 

+19,+21 

0.38 + 0.18 

< 0.12 

< 0.48 

< 1.25 

+21,+24 

0.37 + 0.19 

< 0.11 

< 0.44 

< 1.20 

+24,+27 

0.74 + 0.18 

0.16 + 0.05 

0.64 + 0.20 

0.87 + 0.34 

+40,+46 

7.77 + 0.25 

4.71 + 0.54 

18.82 + 2.15 

2.42 + 0.29 

+47,+50 

5.99 + 0.20 

1.99 + 0.20 

7.95 + 0.79 

1.33 + 0.14 

G29.96-0.02 

+3,+5 

1.55 + 0.59 

0.53 + 0.17 

2.12 + 0.70 

1.37 + 0.69 

+7,+ll 

20.36 + 6.13 

16.31 + 0.47 

65.22 + 1.86 

3.20 + 0.97 

+11,+15 

6.06 + 1.85 

2.22 + 0.23 

8.88 + 0.92 

1.47 + 0.32 

+16,+18 

1.43 + 0.67 

0.37 + 0.14 

1.48 + 0.56 

1.03 + 0.60 
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Table 3—Continued 


Clsr Range 

(km s~^) 

iV(HF) 

(10^^ cm“^) 

iV(p-H20) 

(10^^ cm“^) 

iV(H20)tot 

(10^^ cm“2) 

Af(H20)tot 

Ar(HF) 

G29.96-0.02 

+51,+56 

3.79 ± 0.92 

0.56 ± 0.13 

2.22 ± 0.52 

0.59 ± 0.20 

+57,+62 

8.18 ± 2.30 

2.35 ± 0.35 

9.41 ± 1.40 

1.15 ± 0.37 

+65,+73 

22.49+ 2.44 

9.17 ± 0.59 

36.70 ± 2.35 

1.63 ± 0.21 

+74,+77 

1.80 ± 0.90 

0.22 ± 0.08 

0.90 ± 0.30 

0.50 ± 0.30 

+89,+93 

2.75 ± 0.84 

1.67 ± 0.55 

6.67 ± 2.18 

2.42 ± 1.08 

G34.3+0.1 

+8,+10 

0.57 ± 0.19 

< 0.25 

<1.00 

< 1.74 

+10,+13 

15.19 ± 5.17 

9.66 ± 0.46 

38.62 ± 1.85 

2.54 ± 0.87 

+13,+16 

3.13 ± 0.16 

1.11 ± 0.16 

4.43 ± 0.62 

1.42 ± 0.21 

+17,+26 

2.36 ± 0.53 

0.43 ± 0.32 

1.71 ± 1.28 

0.73 ± 0.57 

+27,+30* 

12.63 ± 0.39 

7.79 ± 1.22 

31.17 ± 4.87 

2.47 ± 0.39 

+40,+44 

2.00 ± 0.19 

2.34 ± 0.47 

9.36 ± 1.89 

4.68 ± 1.04 

+44,+47 

6.54 ± 0.35 

4.33 ± 0.64 

17.32 ± 2.56 

2.65 ± 0.42 

+48,+50 

13.26 ± 0.33 

4.05 ± 0.52 

16.20 ± 2.08 

1.22 ± 0.16 

W33A 

+20,+23 

0.51 ± 0.09 

0.29 ± 0.18 

1.15 ± 0.71 

2.27 ± 1.46 

+24,+27 

7.05 ± 0.38 

3.61 ± 0.51 

14.45 ± 2.05 

2.05 ± 0.31 

+43,+46 

1.24 ± 0.14 

0.92 ± 0.37 

3.69 ± 1.49 

2.97 ± 1.25 

W31C 

+41,+47 

>23.4 

18.58 ± 0.35 

74.32 ± 1.39 

<3.2 
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Table 3—Continued 


Clsr Range 

(km s~^) 

iV(HF) 

(10^^ cm“^) 

iV(p-H20) 

(10^^ cm“^) 

iV(H20)tot 

(10^^ cm“^) 

iV(H20)tot 

Af(HF) 

DR21(OH) 

-11,-9 

0.52 ± 0.12 

<0.12 

< 0.47 

< 0.90 

+11,+13* 

4.36 ± 0.17 

0.90 ± 0.11 

3.59 ± 0.43 

0.83 ± 0.10 

+ 13,+16 

3.43 ± 0.37 

1.03 ± 0.72 

4.11 ± 2.89 

1.20 ± 0.85 

W3(OH) 

-22,-18 

1.70 ± 0.14 

0.19 ± 0.10 

0.76 ± 0.40 

0.45 ± 0.24 

-13,-8 

6.61 ± 0.32 

4.90 ± 0.29 

19.61 ± 1.14 

2.97 ± 0.25 

-6,-4 

1.34 ± 0.22 

0.47 ± 0.10 

1.89 ± 0.38 

1.40 ± 0.37 

-4,-2 

1.61 ± 0.14 

0.77 ± 0.12 

3.08 ± 0.47 

1.91 ± 0.34 

-2,+0 

2.58 ± 0.16 

0.56 ± 0.23 

2.25 ± 0.93 

0.87 ± 0.37 

+0,+2 

8.19 ± 0.11 

6.78 ± 0.08 

27.13 ± 0.31 

3.31 ± 0.06 

W3 IRS 5 

-22,-17 

7.05 ± 0.63 

3.20 ± 0.24 

12.82 ± 0.94 

1.82 ± 0.21 

-6,-3 

1.56 ± 0.34 

0.51 ± 0.23 

2.06 ± 0.92 

1.32 ± 0.66 

-2,+2 

3.61 ± 0.45 

1.04 ± 0.38 

4.17 ± 1.50 

1.15 ± 0.44 

G-0.02-0.07 

-46,-37 

17.34 ± 1.78 

12.34 ± 9.89 

49.37 ± 39.57 

2.85 ± 2.30 

+ 13,+18 

5.17 ± 0.66 

2.14 ± 1.44 

8.55 ± 5.74 

1.65 ± 1.13 

G-0.13-0.08 

-45,-41 

6.01 ± 1.21 

1.52 ± 1.15 

6.09 ± 4.58 

1.01 ± 0.86 

+20,+25 

12.26 ± 1.89 

6.47 ± 2.68 

25.88 ± 10.74 

2.11 ± 0.93 










Table 3—Continued 


Vlsr Range 

iV(HF) 

iV(p-H20) 

iV(H20)tot 

Af(H20)tot 

V(HF) 

(km s~^) 

(10^^ cm“^) 

(10^^ cm“^) 

(10^^ cm“^) 


+26,+30 

4.36 ± 1.56 

1.82 ± 0.90 

7.26 ± 3.62 

1.66 ± 1.02 


®'iV(H 20 )tot was derived from our iV(p-H 20 ) measurements assuming an ortho¬ 
para ratio of 3 for water (Flagey et al. 2013). 

*Our iV(p-H20) value is higher by up to a factor 3 than that reported by Flagey 
et al. (2013) toward W49N and G34.3-1-0.1; our iV(p-H20) value is lower by ~ 
50% than that reported by Flagey et al. (2013) toward DR21(OH). Details can 
be found in Section 6.2. 








